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ABSTRACT 

We investigate whether a circumbinary gas disc can coalesce a supermassive black hole 
binary system in the centre of a galaxy. This is known to be problematic for a prograde disc. 
We show that in contrast, interaction with a retrograde circumbinary disc is considerably 
more effective in shrinking the binary because there are no orbital resonances. The binary 
directly absorbs negative angular momentum from the circumbinary disc by capturing gas 
into a disc around the secondary black hole, or discs around both holes if the binary mass ratio 
is close to unity. In many cases the binary orbit becomes eccentric, shortening the pericentre 
distance as the eccentricity grows. In all cases the binary coalesces once it has absorbed the 
angular momentum of a gas mass comparable to that of the secondary black hole. Importantly, 
this conclusion is unaffected even if the gas inflow rate through the disc is formally super- 
Eddington for either hole. The coalescence timescale is therefore always ~ M2/M, where Mo 
is the secondary black hole mass and M the inflow rate through the circumbinary disc. 
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1 INTRODUCTION 

Astronomers now generally agree that the centre of every reason- 
ably large galaxy contains a supermassive black hole (SMBH). 
Moreover the mass of this hole correl ates (at least at low red- 
shift) with properties of the host galaxy dFerrarese & Merritt 20001 : 
iGebhardt et al. 20001 : iHiiring & Rix 20041) . In the hierarchical pic- 
ture of structure growth, small galaxies merge to produce large 
ones, promoting accretion on to their central SMBHs, and appar- 
ently causing these holes to coalesce. 

The favoured mechanism for driving the holes closer is dy- 
namical friction. However it is unclear that this can bring them 
close enough for gravitational wave losses to complete the co- 
alescence, since the frictional process itself scatters away the 
stars causing it, and refilling of the loss cone is apparently too 
slow. This is often called the 'final parsec problem', as dynami- 
cal friction typically stalls at s uch separations between the holes 
dMilosavlievic & Merritt. 20031) . 

A possible way of overcoming this problem is interac- 
tion with gas orbiting in a di sc just outside the SMBH binary 
( Armitage & Nataraian 20051: MacFadven & Milosavlievic 20081 : 



lLodato et al. 20091 ; ICuadra et al. 2009fT There has also been much 
discussion of cases where an SMBH binary is embedded in a disc 
jEscala et al. 2005l ; lDotti et al. 20071 2009). It is implicitly assumed 
that dissipative torques make the disc coplanar with the binary. 
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Studies of the circumbinary disc problem have so far considered 
prograde discs, i.e. those rotating in the same sense as the bi- 
nary. Then tidal interaction with the binary turns the disc into a 
decretion disc, which transports a ngular momentum outward, but 
with little inward mass transport jLodato et al.. 20091) . If the disc 
mass is large enough to carry away the binary angular momentum, 
a decretion disc is vulnerable to the self-gravitational instability 
( lLodato et al. . 20091) . This can rob the disc of the gas it needs to 
drive further angular momentum loss, halting the binary shrinkage. 

This makes it doubtful that a prograde disc can ever in 
practice shrink the binary separation from a ~ 1 pc to the 
point (a ~ 10~ 2 pc) where gravitational wave losses can drive 
it to coalescence. However the separation of the SMBH bi- 
nary is much smaller than the interacting galaxies themselves 
so it is highly unlikely that the central gas flows are al- 
ways prograde. These flows also receive randomly-directed in- 
jections of energy and momentum from star formation and su- 
pernovae, suggesting that retrograde flows are as likely as pro- 
grade. This kind of chaotic accretion gives a plausible picture 
of the mass and spin evolution of a central accreting b lack hole 
dKing & Pringle 2006|lKing et al. 2008llHobbs et al. 2010h (but see 
iBerti & Volonteri 20081 ; iMaver et al. 20071 for an alternative pic- 
ture). It also implies that if nothing else manages to drive an SMBH 
binary to coalescence, it is highly likely that at some point there will 
be a retrograde coplanar disc surrounding the binary. This situation 
would arise if for example an earlier episode of prograde gas accre- 
tion failed to coalesce the binary, and was followed by a retrograde 
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accretion event. As we will see, retrograde discs behave quite dif- 
ferently from prograde ones, and may offer a solution to the final 
parsec problem. 



2 PROGRADE VERSUS RETROGRADE 

We start by contrasting the main qualitative features of the prograde 
and retrograde cases. These stem from the physics of the interac- 
tion between the binary and the disc, where dissipative torques try 
to share the angular momenta. If the disc is prograde, this interac- 
tion shrinks the binary but moves the inner edge of the disc out- 
wards, reducing the torque shrinking the binary. If instead the disc 
is retrograde, the effect is to shrink the binary, but also to move the 
inner edge of the disc inwards (see Appendix [At. 

A prograde disc becomes a decretion disc, transporting 
angular momentum outwards with little net mass transport 
dLodato et al.. 20091) . A retrograde disc instead remains an accre- 
tion disc, transporting angular momentum outwards and mass in- 
wards. As we shall see, the long-term evolution of the disc-binary 
system is radically different in the two cases. 

In particular, the disc-binary torque is quite different. For 
a prograde disc the tidal interaction occurs mainly through reso- 
nances. These occur when 



SI 2 = m 2 (Q. - cj) 2 , 



(1) 



where <±> is the binary orbital frequency, Q(R) is the disc angular 
velocity and m = 1, 2, ... is the wave mode n umber (for example 
see the analysis in lPapaloizou & Pringle 19771) . So there are reso- 
nances at radii where 



m+l 



(2) 



where £1 and u> have the same sign. Resonances outside the binary 
orbit (i.e. with Q < oi) correspond to the positive sign in the de- 
nominator of equation [2] and so appear at radii where 



Q.(R) 



1 2 3 

2' 3' 4' 



(3) 



The dominant interaction then involves the 2:1 (more strictly 1:2) 
resonance. 

By contrast, in a retrograde disc CI and to have opposite signs, 
and equation Q] requires |Q| > \co\, so there are no resonances in 
a circumbinary retrograde disc. The disc-binary interaction is di- 
rect, as the inner edge of the disc starts to impinge on the sec- 
ondary black hole. A retrograde circumbinary disc remains an ac- 
cretion disc whose material is gravitationally captured by the bi- 
nary, directly reducing its angular momentum. This is inherently 
more promising for shrinking the binary towards coalescence than 
the prograde case, where the binary dams up the disc. 

It is important to understand that 'capture' simply means that 
the gas orbits a particular hole, and so has added its (negative) an- 
gular momentum to the binary orbit. It does not imply that the rel- 
evant hole must actually accrete this gas (although it may). Once 
captured the gas is bound to the hole and thus may be treated as a 
single body. Some or all of this captured gas can be expelled, for ex- 
ample by radiation pressure. Provided that this process is isotropic 
in the frame of the hole this does not change its orbital angular mo- 
mentum, and so has very little effect on the orbital dynamics and 
eventual coalescence (cf equation[22l below). 



3 WHERE DOES THE MASS GO? 

The effect of a retrograde circumbinary disc differs in detail de- 
pending on how the captured mass is distributed between the two 
black holes. Accordingly we look at the reaction of test particles to 
the binary. To make things simple we first consider a circular binary 
with a low mass ratio, i.e. MijM\ < 0. 1. Then to first order we can 
treat the primary as fixed and the secondary as following a circular 
orbit of radius a around it with velocity V = (GM]/a)'' 2 . The cir- 
cumbinary disc gas is on circular orbits with velocity =s (GM] /R) l/2 
at each radius R. We assume that the disc creeps slowly inwards 
from a large radius by the usual viscous evolution. As it is counter- 
rotating we can ignore all resonant effects and assume everything 
is ballistic until orbits begin to cross and fluid effects appear. We 
define the effective radii Ri,R2 of the two holes as the radii within 
which gas particles are captured by each hole, e.g. by forming a 
disc around one or other of them, so that this captured gas has 
the same net specific angular momentum as the relevant hole. Evi- 
dently Ri,R2 cannot in practice be larger than the individual Roche 
lobes for each hole. 

For the disc particles orbiting closest to M 2 the interaction 
with M 2 is initially hyperbolic and so we can use the impulse ap- 
proximation. Here the relative velocity is approximately 2V. If the 
disc edge is at a radius R = a + b, where b <k a, the impulse ap- 
proximation shows that the disc particle acquires an inward radial 
velocity 



(GM 2 \ 2b 



2b GM 2 



bV 



(4) 



The inner edge of the circumbinary disc is not significantly per- 
turbed if its distance b from the secondary's orbit is large enough 
that U R < V, i.e. V 2 > GM 2 /b, or equivalently b > {M 2 IM{)a. 
From this we conclude that the secondary cleanly pulls the gas from 
the unperturbed inner edge of the disc at b provided that R 2 ^ b, 
i.e. 



R 2 > M 2 
a ~ M\ 

Comparing this with the Roche lobe constraint 

^0.4/(^ 
a J M| 



(5) 



(6) 



where / < 1 is a dimensionless factor, we see that the secondary 
captures almost all the gas for mass ratios 



M 2 
~M~, 



< 4 crit = 0.25/ 3/2 = 0.25. 



(7) 



Equivalently we can define a Safronov number dSafronov. 19721) as 
used in discussing accretion of planetesimals from a protoplanetary 
disc: 







M 2 a 
~M~\~R 2 



(8) 



where u csc ,Uorb are the escape velocity from the secondary's effec- 
tive radius, and its orbital velocity respectively. This measures how 
much the gas is gravitationally perturbed before being captured. 
Small © implies little perturbation. We see that the criterion for 
secondary capture is just < 1, that is, the gas is cleanly captured 
without significant perturbation. 

For larger mass ratios the flow becomes more complex and 
it is likely that some gas falls towards the binary centre of mass, 
producing some form of primary capture. In numerical simulations 
(Section we will find that even in this case most of the gas is 
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captured by the secondary. So in general the secondary captures 
most of the mass, and only in a major merger with q > g crit does 
the primary capture a non-negligible amount from a retrograde cir- 
cumbinary disc. 



4 ECCENTRICITY GROWTH 

A retrograde circumbinary disc can decrease both the energy and 
angular momentum of the SMBH binary, and so change its eccen- 
tricity. A simple argument shows how this happens. We consider 
a slightly eccentric binary orbit (the orbit is never exactly circu- 
lar if it is shrinking) and again for simplicity assume that the mass 
ratio is sufficiently extreme that we can regard the primary black 
hole as effectively fixed, and only the secondary as interacting with 
the disc. None of these restrictions will affect our conclusions (see 
Appendix [A] for a fuller discussion). 

Momentum conservation shows that capture of disc gas al- 
ways reduces the secondary's orbital velocity (see equation |18t . 
This holds both for a direct collision, or (more commonly) if the 
secondary captures gas into a bound disc around itself. The sec- 
ondary's mass cannot decrease in these interactions. Its specific 
orbital energy therefore always drops, so that the binary semima- 
jor axis a decreases. But for mass capture into a symmetrical disc 
around the secondary near apocentre, the new orbit must retain the 
same apocentre as the old one (it must pass through this point, and 
the radial velocity remains zero there). Given the decrease of the 
semimajor axis, this means that capture near apocentre tends to in- 
crease the orbital eccentricity e (since a[l + e] remains constant). 
The eccentricity is evidently 

AM 

e — — (9) 
M 2 

where AM is the amount of mass the secondary has captured. Ex- 
actly the same reasoning shows that for mass capture near peri- 
centre, the quantity a(l — e) has to stay constant despite a further 
decrease in a - in other words, the eccentricity must decrease here 
by about the same amount (AM/Mj) it increased at apocentre. 

The secondary obviously captures at all points in between 
apo- and pericentre, but the effects are opposed for significant 
times. If the eccentricity is initially small these times are nearly 
equal, and e stays small as the orbit shrinks, provided that AM/M 2 
is below a certain value we will derive shortly. 

If on the other hand the orbit is initially quite eccentric, or 
the mass grows significantly in one orbit, the pericentre may be 
too small to allow disc interaction. Then e grows with a{\ + e) = 
Oq - constant, and the pericentre distance goes as 



p = a{\ - e) = 2a - a . 



(10) 



The binary therefore coalesces once the original semimajor axis 
has halved. We show below (equation [22} that this occurs once the 
secondary has absorbed the (negative) angular momentum of disc 
gas with mass comparable to its own. If gas flows inwards through 
the circumbinary disc at the rate M the timescale for coalescence is 



Mi 
M 



(11) 



We stress again that this does not require either hole to accrete this 
gas but only to capture the gas into a bound orbit around the hole. 
In particular, M can formally exceed the Eddington accretion rate 
for either hole without affecting the orbital shrinkage. If the accre- 
tion rates on to either hole were super-Eddington then the captured 
gas would be blown away by radiation pressure from the disc(s) 



around the secondary (and possibly primary) black hole, without 
significantly changing the binary orbital evolution (cf equation[22l. 

The critical eccentricity separating cases where the binary re- 
mains almost circular from those of growing eccentricity depends 
on the surface density distribution of the circumbinary disc. In a 
real three-dimensional disc with scale height H and aspect ra- 
tio H/R the surface density tails off over a length-scale H(a) ~ 
(H/R)a. This suggests that the critical eccentricity dividing these 
two cases is just 



H 



(12) 



Thus any binary starting with e > e crit , or achieving it by capturing 
a large mass (comparable to the secondary's) in one orbit, must be- 
come very eccentric. We note that even a preceding episode of pro- 
grade ac cretion can leave the binary with an eccentricity exceeding 
e crit (e.g. lCuadra et al. 200 9;!), so that growth to high eccentricity is 
very likely if accretion is chaotic. 

Our conclusion s about eccentricity growth agree with those 
of iDotti et al. (20091) . who considered a related but different prob- 
lem. A secondary black hole was injected with significant eccen- 
tricity into a pre-existing dense circumnuclear disc surrounding a 
primary black hole. The secondary's orbit was initially retrograde 
with respect to this interior disc. However the secondary was able 
to interact with enough gas in less than one orbit that it cancelled all 
of its angular momentum. The secondary then briefly had approxi- 
mately zero angular momentum before capturing more gas and so 
changing its angular momentum to prograde. In this paper we re- 
strict ourselves to accretion events on much smaller scales than in 
IDotti et al. (20091) . This difference in lengthscales is important. At 
the smaller scales we consider, a disc with mass Mj » Mi would 
probably be self-gravitating, and this change of angular momentum 
sign is unlikely to occur. Our simulations agree with this conclu- 
sion: once interaction with the disc has cancelled the orbital angular 
momentum (and thus much of the orbital energy) of the secondary, 
the binary coalesces. 



5 ORBITAL EVOLUTION WITH A RETROGRADE 
CIRCUMBINARY DISC 

We can now make analytic estimates of the orbital evolution as 
the binary interacts with an exterior disc. For simplicity we again 
assume that q = MijM\ « 1, so that the secondary has specific 
orbital energy and angular momentum E, J 2 where 



E = 
and 



GM, 
la 



1 



-MrY 



GM, 



r = GMia{\ - e-) 



(13) 



(14) 



We have seen above that in many cases the binary eccentricity 
grows quite strongly. In the limit the secondary interacts with the 
disc only very near apocentre r = a{\ + e). Here its velocity v ap is 
purely azimuthal, with 



[a(l + e)v, f f = J- 

which by equation[74l gives 

2 _ GMi 1 - e 
~ ~a~TT~e 



(15) 



(16) 



Near apocentre the secondary interacts with disc material moving 
with azimuthal velocity va x < 0, with 



4 C.J. Nixon, P.J. Cossins, A.R. King & J.E. Pringle 



GMi 



GM, 



r a(l + e) (17) 
We assume that a mass AM of disc matter is captured into 
orbit about the secondary near apocentre, as we discussed above, 
so that all of its orbital angular momentum is transferred to the 
secondary. To allow for mass loss from the subsequent accretion 
process on to this hole (if for example this is super-Eddington, or 
mass interacts gravitationally with the secondary but is all accreted 
by the primary) we assume that the effective mass of the hole plus 
disc becomes M 2 + a AM, with < a < 1. Then conservation of 
linear momentum gives 



M 2 u a p - AMv d 



(M 2 + aAM)u 



(18) 



where u is the new apocentre velocity of the secondary plus its 
captured gas disc. The changes AE, Aa in orbital specific energy 
and semi-major axis are given by 



GM\ 1 , 1 2 

-Afl = At = —U v~ 

2a 2 2 2 Jp 



Combining eq uations 1 1 61 1 7 1 1 8l and 1 1 9 1 gives 
Aa -2 AM , 



a(l + e) M ; 



■[(l-e) 1/2 + a(l-e)] 



(19) 



(20) 



to lowest order in AM. We noted above that if the interaction is 
confined to the immediate vicinity of apocentre then the apocentre 
distance cio = a{\ + e) stays constant in the subsequent evolution. 
Then 1 + e = a /a and 1 - e = 2 - a /a, so that the LHS of equa- 
tion[20]is simply A(l - e)/ao and therefore equation [20]becomes 



(21) 



2AM , „ 

A(l-e) = - — -[(1 - ef' 2 + a(l - e)]. 

jVl2 

Using M2 = M20 + aM, where M20 was the mass of the sec- 
ondary hole when e was zero, and M is the total mass since trans- 
ferred from the disc, this integrates to give 

-1 , (22) 



v M 2 o + aM 

Hence in this approximation the binary coalesces (i.e. 1 - e = 0) 
once the disc has transferred a mass equal to the secondary black 
hole, i.e. after a time f co (cf equationUTJ, independently of the frac- 
tion a. This means that mass loss has no effect in slowing the in- 
spiral. We note that equation [20] implies that the energy dissipated 
in the shrinkage of an eccentric binary is 

GMi AM 



- M 2 AE 



2a 



(23) 



per binary orbit, which is less than that produced by viscous dissi- 
pation in the disc (this has to pull in a mass larger than AM on each 
orbit). 

All other cases give similar timescales ~ f co for coalescence. 
For example if the orbit stays circular we can set e = in equa- 
tionl20land find 
Aa AM 

— = -2(1 + a) (24) 
a M 2 

so that a oc M~ 2(i+a \ If the secondary gains the transferred mass 
we have a = 1, a oc M^ 4 , while if the transferred mass (but not the 
angular momentum) ends up on the primary we have a = 0, a oc 
Mj 2 . Shrinking the binary from a = 1 pc to a = 10~ 2 pc, where 
gravitational radiation rapidly coalesces it, requires the transfer of 
between 2 and 9 times the original mass of the secondary in the two 
cases. This is larger than in the eccentric case because the torque on 
the binary now decreases with a (indeed the required mass would 
be formally infinite without gravitational radiation). 



6 SIMULATIONS 

We now use Smoothed Particle Hydrodynamics (SPH) simulations 
to compare with the analytical arguments above. 



6.1 Code setup 

We use a fully 3D conservative L agrangian implementation o f 
the SPH al g orithms, for ex ample ISpringel & Hernquist (20021) ; 
Price (20051) : iRosswog (20091) . We neglect gas self-gravity and so 
are able to use link ed-lists of part icles rather than the usual tree for 
neighbour finding JPeegan. 200^) . We assume an isothermal disc, 
i.e. P = cjp, where the temperature and hence sound speed c s is 
constant for all particles at all times. We make this choice for sim- 
plicity, and anticipate that the interaction of gas with the binary is 
not greatly affected by it. In particular we noted in Section[5]above 
(cf equationl23l that binary shrinkage does not greatly increase the 
heating of the disc. We integrate a 'live' binary which feels the 
back reaction of the gas and generates self-consistent orbits for the 
binary and the gas. 

Our aim here is to understand the dynamical interaction of 
the disc with the binary, rather than long-term behaviour such 
as the timescales for coalescence, which we know depends on 
the total mass transferred through the disc (see Section [5] above). 
Accordingly we use only the standard SPH artificial viscosity 
jRosswog. 200 9) and do not include a physical viscosity. Our simu- 
lations run only for a few tens of binary orbits, so that the viscosity 
scheme employed should not affect our results. 

Our code units are dimensionless and take the initial separa- 
tion and mass of the binary as unity. The period of the binary in 
its initial configuration is thus 2n. The general disc setup is gas, 
of mass M d , spread from R m to R ollt with a surface density follow- 
ing a power law in radius, i.e. S = £0 (R/Rq) p with p typically 
= 1. The initial vertical structure in the disc is a Gaussian. We set 
the sound speed in the gas to a value ensuring that our neglect of 
self- gravity is justifi ed. Thus we arrange that the Toomre Q param- 
eter dToomre. 1964) exceeds a minimum value (> 1) throughout 
the disc. Typically we arrange Q > 5 for a disc of mass M rf = 10~ 2 . 
The particles are all on initially circular orbits. 

The simulations detailed below were repeated with half a mil- 
lion, one million and two million particles and were deemed to have 
converged. All of the simulations detailed below initially contain 
one million particles. 1 



6.2 Where does the mass go? 

Two main parameters govern the gas flow for a retrograde disc- 
binary system: the binary mass ratio and the 'capture' radius of 
each binary component. As we discussed in Section [3] with mass 
ratios g - 1 , both the primary and secondary can interact with the 
inner edge of the disc so we should expect gas capture by both 
objects. However if q <K 1, an inward-moving disc reaches the 
secondary before the primary. Then the secondary's effective ra- 
dius may be large enough for it to capture all of the gas. If instead 
this radius is small the gas may be perturbed towards the primary 
rather than captured. In Section [3] we estimated the critical value 
dividing these cases (equation[5)- We study the effects of changing 
parameters here. 



1 Movies of these simulations are available and can be found at: 
www . astro . le . ac . uk/users/c jn!2/RetroBinaryMovies . html 
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It is of course currently impossible to follow gas accretion 
down to the innermost stable circular orbits of the two black holes 
while simultaneously following the gas flow from the circumbi- 
nary disc. The analytic work above suggests that the secondary 
hole is always surrounded by a gas disc denser than the circumbi- 
nary disc and any gas flow from it, which therefore captures any 
gas impacting it. If initially there had been no gas around the sec- 
ondary, gas entering i? cllplure = aMil2M\ would be captured (but 
not necessarily accreted). This gas would spread to form a disc 
around the secondary, on a timescale short with respect to the 
mass transfer timescale in the circumbinary disc. The secondary's 
disc then forms an obstacle for further inflow from the circumbi- 
nary disc. As further mass is captured in this way, this disc can 
become at most as large as some fraction of the Roche lobe (cf 
equation [6). In the language of planetary dynamics this condition 
is formally equivalent to taking the secondary's Safronov num- 
ber, defined as (planetary escape velocity /orbital velocity) 2 /2, (e.g. 



'tV/< 

I Safronov 19721 ; Hansen & Barman 20071 ) to be of order 2.5q 



From this analysis we expect the capture radius of the holes 
to be of order or smaller than the Roche lobe. For completeness 
we also consider the possibility of still smaller capture radii. Once 
a particle moves within the capture radius it is assumed to have 
impacted upon the disc assumed to be present inside the capture 
radius. This gas could be accreted by the black hole or expelled by 
radiation pressure from the innermost parts of the disc. In either 
case it is clear that once it has been 'captured' it plays no signif- 
icant further role in the binary dynamics. Accordingly we add its 
mass and momentum to the relevant hole and remove it from the 
simulation. 



6.2.1 Secondary capture 

The first simulation has a mass ratio of q = 0. 1 and effective radii 
0.15, 0.2 for the secondary and primary. Equation[5]implies that the 
secondary should capture the inner edge of the circumbinary disc 
without significantly perturbing the remaining disc particles. From 
equation[8]the secondary's Safronov number is 2/3. 

We start with the retrograde disc extending from 1 to 2 in ra- 
dius, and a circular binary. At first the disc inner edge is slightly 
perturbed, but not enough to provide any noticeable accretion onto 
the primary. After one binary orbit the system settles and the sec- 
ondary smoothly captures the inner edge of the disc. Figure[Tjshows 
the state of the simulation at t = 50. At this time, of the particles 
accreted, the secondary has captured ~ 100% (149759 particles) 
and the primary has captured ~ 0% (3 particles). 




Figure 1. Image of the simulation from Section |6.2.1| at time t = 50. The 
binary is represented by the two dots. The axes are in code units with the 
log of the column density given by the bar. 




Figure 2. Image of the simulation from Section l6.2.2l at time t = 50. 



6.2.2 Primary capture 

The second simulation also has a mass ratio of q = 0.1, but this 
time we use an effective radius of 0.05 for the secondary with again 
0.2 for the primary capture radius, so that its Safronov number is 
= 2. Equation ((5} implies that the secondary should perturb the 
disc particles significantly without directly accreting all of them. 

The initial disc is exactly the same as for the simulation in 
Section |6.2.1| The binary is again initially circular. 

At first, the disc inner edge is significantly perturbed, with 
some particles on orbits passing close to the primary. In the first 3 
- 4 orbits the gas flow is very chaotic. After this the system settles 
into a quasi-steady state. Both holes capture from the retrograde 
circumprimary disc (as shown in figure[2]l. The secondary both dis- 
turbs and captures particles from the outer edge of the circumpri- 
mary disc, and perturbs more particles into it. The circumprimary 



disc is warped, eccentric and precessing. This is probably a result 
of particle noise destroying the symmetry about the orbital plane, 
as might indeed happen in a realistic situation. 

At time t = 50 the secondary has captured -95% (147642) 
particles and the primary has captured ~ 5% (8193 particles). So 
even in this case the secondary still takes most of gas captured from 
the circumbinary disc. The binary probably coalesces once the disc 
has transferred a mass ~ M 2 , so if the disc is more massive than 
this, the coalesced hole eventually accretes the remainder. Clearly if 
the secondary's mass or its capture radius is made arbitrarily small, 
it would accrete very little, and most of the mass would be cap- 
tured by the primary. However it is clear that this requires extreme 
choices of these parameters. 
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Figure 3. Images of the simulation from Section l6.2.3l at time ; = 25. The 
upper panel shows the disc face on. The lower panel shows the same plot 
but viewed at an angle of 15° to the plane of the disc. 



6.2.3 Dual Capture 

Here we look at the interaction of an near equal-mass binary with 
the disc. We use q = 0.5 and the same disc as in Sections [6.2.1l and 
16.2.21 Again the binary is initially circular. Here we use capture 
radii of 0.01 for both holes to allow formation of circumprimary 
and circumsecondary discs. 

During the initial 3-4 binary orbits the flow is very chaotic, 
with mass captured by both holes. As more gas is captured (but this 
time not 'accreted') by both holes, circumprimary and circumsec- 
ondary discs are formed which persist throughout the simulation. 
The discs are supplied by streams from the circumbinary disc. We 
show the simulation at time / = 25 in figure refboth. At this time the 
primary has captured 48% (57969) and the secondary 52% (62189) 
of the accreted particles. Together with the previous simulation, this 
shows that the primary can gain significant mass only if the mass 
ratio is close to unity. We note again that the discs are not planar, 
and show a significant tilt with respect to the binary plane. Again 
this is probably a result of particle noise removing the symmetry 
about the orbital axis. Globally angular momentum is conserved 
but locally the streams that supply these discs need not be planar. 
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Figure 4. Eccentricity growth and decay of the simulations in Section 
Time is in code units. The solid line is the first simulation (capture at apoc- 
entre but not at pericentre), the dashed line is the second simulation (capture 
at both pericentre and apocentre) and the dotted line is the third simulation 
(capture at pericentre and not apocentre) 



6.3 Eccentricity growth 

The analytic arguments of Section [4] suggest that capture from a 
retrograde circumbinary disc at apocentre and pericentre of the bi- 
nary orbit increase and decrease its eccentricity respectively. Here 
we show three SPH experiments exploring this. 

The first simulation has an initially circular binary, with q = 
10~ 3 , where the secondary is embedded just inside the inner edge 
of a retrograde disc. In code units the disc is spread radially from 
0.8 to 1.5. The disc has an initial mass Mj = 10~ 2 . 

The second simulation also starts with the secondary embed- 
ded in the disc. However this time the inner edge of the disc extends 
much further in, to 0. 1 . 

The third simulation starts with a retrograde disc interior to the 
binary, i.e. a circumprimary disc extending from 0.1 to 1.0. In this 
case we start with initial eccentricity 0.5 and the binary at apocen- 
tre - so the secondary begins outside the circumprimary disc but 
plunges into it before it reaches pericentre. 

All of these simulations have capture radii 0.1, 0.01 for the 
primary and secondary respectively. 

We show in figure|4]how the eccentricity of the binary evolves 
in the three cases. In the first simulation the binary captures enough 
gas at apocentre on the first orbit to plunge inside the inner edge of 
the disc. This means that the eccentricity growth at apocentre can- 
not be moderated by decay at pericentre. The eccentricity therefore 
grows as shown in figure|4] In the second simulation the secondary 
is still capturing when it reaches pericentre and so the eccentric- 
ity decays. It is notable that the eccentricity never returns precisely 
back to zero in this case. This happens because the secondary cap- 
tures unequal amounts of gas at apo- and pericentre. This is reason- 
able, as its velocity is higher at pericentre and so it has less time to 
interact with the gas, In addition the mass and angular momentum 
at each radius is not constant. 

In the third simulation the secondary can initially only capture 
from the circumprimary disc at pericentre and so the eccentricity 
decays. After approximately 3-4 orbits the secondary has damped 
out most of its initial eccentricity to an almost circular orbit. After 
a short interaction with the inner parts of the disc it reaches its 



inner edge. Here there is a brief interval of eccentricity growth as it 
captures from this inner edge of a (now) circumbinary disc. 

All of the results of these simulations agree with the analytic 
arguments of Section|4] 



7 DISCUSSION 

We have seen that accretion from a retrograde circumbinary disc 
can be considerably more effective in shrinking a supermassive 
black hole binary system than accretion from a corresponding pro- 
grade disc. Coalescence occurs if the secondary black hole captures 
angular momentum from a gas mass comparable to its own. 

The reason for this effectiveness is the absence of or- 
bital resonances in a retrograde disc. This cannot react to 
tidal torques from the binary and become a decretion disc, 
which is what tends to slow the evolution in the prograde 
case llLodato et al. 20091; ICuadra et al. 20091 : iDotti et al. 20071 ; 
iMacFadven & Milosavljevic 20081) . Instead, the disc directly feeds 
negative orbital angular momentum to the secondary black hole. 
An important aspect here is that there is no restriction on the rate 
at which this can occur. In particular gas can flow inwards at rates 
higher than Eddington for even the primary black hole without any 
significant effect on the amount of gas required to coalesce the bi- 
nary. 

Capture to the primary hole is negligible unless the mass ratio 
is very close to unity, and is never dominant (even if 'capture' ac- 
tually leads to 'accretion', which is not required). This may be im- 
portant in interpreting attempts to observe merger events. In many 
cases the binary may develop very high eccentricity before gravi- 
tational wave emission coalesces it. The latter rapidly damps the 
eccentricity during the final inspiral, but leaves a residual value 
which is likely to be significant in any merger event detectable by 
gravitational wave observatories such as LISA. If the secondary 
hole actually accretes the captured mass it will acquire a large spin 
(Kerr parameter a 1) antiparallel to its binary orbit, which would 
be potentially detectable in the LISA waveform. Moreover if the 
mass ratio is very close to unity it is conceivable that the primary 
hole might do the same. This could in principle favour high black 
hole spin as a result of major mergers (e.g. in giant ellipticals), as 
sometimes proposed. However we should recall that the holes may 
well not gain much mass, particularly if this is captured at super- 
Eddington rates, and the primary definitely does not gain mass (and 
thus spin up) unless the mass ratio is very close to unity. For these 
reasons it seems unlikely that coalescences induced by retrograde 
gas flows produce rapidly spinning merged holes except in rare 
cases. 

We caution finally that so far we have not considered two 
important effects. First, we assumed that the retrograde disc was 
coplanar with the binary, whereas in reality counter-alignment 
must occur over a viscous dissipation timescale. Second, self- 
gravity can deplete the cir cumbinary disc of g as and reduce its abil- 
ity to shrink the binary (cf. lLodato et al. 20091 in the prograde case). 
Since coalescence by a retrograde disc requires M d > M 2 , and self- 
gravity effects appear unless M d < (H/R)(Mi + M 2 ), where H/R is 
the disc aspect ratio, it appears that the last parsec problem is so far 
alleviated only for mergers with M 2 /M\ < H/R. This paper does 
suggest how coalescence might work for larger mass ratios. The 
secondary hole has to absorb negative angular momentum from a 
gas mass comparable to its own, preferably in an eccentric orbit. 
Although gas self-gravity is inevitably important in such an event, 
a retrograde flow has the advantage that there is no limit on the 
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mass inflow rate. We shall investigate these two effects in future 
work. 
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APPENDIX A: BINARY-DISC INTERACTION 

To understand the interaction between the gas and the binary we 
consider the case where the disc particles can be perturbed before 
impacting upon the secondary. The relative velocity of a disc par- 
ticle before interacting is 2V. After the interaction it gains a radial 
velocity (equation |4). In the frame of M 2 the energy of the particle 
is conserved and so 

(2V) 2 = U 2 R + Uj. (Al) 
so that 



8 C.J. Nixon, P.J. Cossins, A.R. King & J.E. Pringle 



Ur * 2V 



G 2 M, 



2 

1 2 - 

b'-V* 



(A2) 



Back in the inertial frame, this implies that after the interaction the 
radial velocity of the disc particle is U R and the azimuthal velocity 
is Up = U T - U, i.e. 

G 2 M 2 

U * = V --b-^T- CA3) 

The particle was on an orbit with eccentricity e = 0, specific energy 
E = GM[/a and specific angular momentum h = (GM2R) 112 . After 
the interaction the particle's specific kinetic energy is T where 

, Uj U\ V 2 2G 2 M\ 

T' = — + — = — + (A4) 

2 2 2 b 2 V 2 { ' 

Thus after the interaction the specific energy of the particle orbit is 
increased to 

E' = E + AE, (A5) 
where 

G 2 M 2 

AE = ^vr (A6) 

and the particle's specific angular momentum is decreased to 

h' =h- Ah, (A7) 

where 

G 2 M 2 

M = a- IV ±. (A8) 

We note that the angular momentum of the particle has the 
opposite sign to that of the secondary, as it must. 

The perturbed particle now has semi-major axis larger than a, 
by an amount Aa = a (AE/E). Thus 

G 2 M 2 

Aa = a-^1, (A9) 
and also for comparison 

Aa a 3 Ml 
b b 1 Mf 

The eccentricity of the particle's new orbit is given by 

Ah AE Ah Aa 

ex = — + — . (All) 

h 2E h 2a 

So we find 

G 2 Ml 

From this we conclude that the interaction with the binary in- 
creases the particle's orbital energy and decreases its orbital angular 
momentum, leading to an eccentric orbit. Now there are two pos- 
sibilities for the disc. Either the disc now becomes eccentric itself. 
Or the orbits of the perturbed particles begin to cross and therefore 
the particles collide, lose energy and share angular momentum. 

As the disc can radiate away the excess energy associated with 
the eccentric disc orbits (cf equation [23} it is clear that the overall 
effect on the disc of the interaction with the binary is to shrink the 
disc. Similarly for the binary we conclude that the binary loses both 
orbital energy and angular momentum and hence shrinks. 



